Black Strings, Low Viscosity Fluids, and Violation of Cosmic Censorship 
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We describe the beliavior of 5-dimensional black strings, subject to the Gregory-Laflamme insta- 
bihty. Beyond the hnear level, the evolving strings exhibit a rich dynamics, where at intermediate 
stages the horizon can be described as a sequence of 3-dimensional spherical black holes joined by 
black string segments. These segments are themselves subject to a Gregory-Laflamme instability, 
resulting in a self-similar cascade, where ever-smaller satellite black holes form connected by ever- 
thinner string segments. This behavior is akin to satellite formation in low-viscosity fluid streams 
subject to the Rayleigh-Plateau instability. The simulation results imply that the string segments 
will reach zero radius in finite asymptotic time, whence the classical space-time terminates in a 
naked singularity. Since no fine-tuning is required to excite the instability, this constitutes a generic 
violation of cosmic censorship. 

PACS numbers: 04.50. Gh, 04.20. q, 04.25.D 



Introduction: While stationary black holes in 4 space- 
time dimensions (4D) are stable to perturbations, higher 
dimensional analogues are not. Indeed, as first illustrated 
by Gregory and Laflamme in the early 90s black 
strings and p-branes are linearly unstable to long wave- 
length perturbations in 5 and higher dimensions. Since 
then, a number of interesting black objects in higher di- 
mensional gravity have been discovered, many of them 
exhibiting similar instabilities (see e.g. [2|). 

An open question for all unstable black objects is what 
the end-state of the perturbed system is. For black 
strings [l[ conjectured that the instability would cause 
the horizon to pinch-off at periodic intervals, giving rise 
to a sequence of black holes. One reason for this con- 
jecture comes from entropic considerations: for a given 
mass per unit length and periodic spacing above a critical 
wavelength Ac, a sequence of hyperspherical black holes 
has higher entropy than the corresponding black string. 
Classically, event horizons cannot bifurcate without the 
appearance of a naked singularity Thus, reaching 
the conjectured end-state would constitute a violation of 
cosmic censorship, without "unnatural" initial conditions 
or fine-tuning, and be an example of a classical system 
evolving to a regime where quantum gravity is required. 

This conjecture was essentially taken for granted un- 
til several years later when it was proved that the gen- 
erators of the horizon can not pinch-off in finite afhne 
time From this, it was conjectured that a new, non- 
uniform black string end-state would be reached Q ■ Sub- 
sequently, stationary, non-uniform black string solutions 
were found however, they had less entropy than 

the uniform string and so could not be the putative new 
end-state, at least for dimensions lower than 13 0- 

A full numerical investigation studied the system be- 
yond the linear regime [8|, though not far enough to 
elucidate the end-state before the code "crashed". At 



that point the horizon resembled spherical black holes 
connected by black strings, though no definitive trends 
could be extracted, still allowing for both conjectured 
possibilities: (a) a pinch-off in infinite affinc time, (b) 
evolving to a new, non- uniform state. If (a), a question 
arises whether pinch-off happens in infinite asymptotic 
time; if so, any bifurcation would never be seen by out- 
side observers, and cosmic censorship would hold. While 
this rnight be a natural conclusion, it was pointed out 
in [To] that due to the exponentially diverging rate 
between affinc time and a well-behaved asymptotic time, 
pinch-off could occur in finite asymptotic time. 

A further body of (anecdotal) evidence supporting the 
GL conjecture comes from the striking resemblance of 
the equations governing black hole horizons to those de- 
scribing fluid flows, the latter which do exhibit insta- 
bilities that often result in break-up of the fluid. The 
fluid/horizon connection barkens back to the membrane 
paradigm 111, and also in more recently developed cor- 



respondences [lj,ll3l. In [Ij it was shown that the dis- 
persion relation of Rayleigh-Plateau unstable modes in 
hyper-cylindrical fluid flow with tension agreed well with 
those of the GL modes of a black string. Similar behavior 
was found for instabilities of a self-gravitating cylinder of 
fluid in Newtonian gravity [l5l |. In [l6l |. using a pertur- 
bative expansion of the Einstein field equations [isj to 
related the dynamics of the horizon to that of a viscous 
fluid, the GL dispersion relation was derived to good ap- 
proximation, thus going one step further than showing 
analogous behavior between fluids and horizons. 

What is particularly intriguing about fluid analogies, 
and what they might imply about the black string case, 
is that break-up of an unstable flow is preceded by for- 
mation of spheres separated by thin necks. For high vis- 
cosity liquids, a single neck forms before break-up. For 
lower viscosity fluids, smaller "satellite" spheres can form 
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in the necks, with more generations forming the lower 
the viscosity (see [l^l for a review). In the membrane 
paradigm, black holes have lower shear viscosity to en- 
tropy ratio than any known fluid [isj . 

Here we revisit the evolution of 5D black strings using 
a new code. This allows us to follow the evolution well 
beyond the earlier study [1]. We find that the dynamics 
of the horizon unfolds as predicted by the low viscosity 
fluid analogues: the string initially evolves to a configu- 
ration resembling a hyperspherical black hole connected 
by thin string segments; the string segments are them- 
selves unstable, and the pattern repeats in a self-similar 
manner to ever smaller scales. Due to finite computa- 
tional resources, we cannot follow the dynamics indefi- 
nitely. If the self-similar cascade continues as suggested 
by the simulations, arbitrarily small length scales, and in 
consequence arbitrarily large curvatures will be revealed 
outside the horizon in finite asymptotic time. 
Numerical approach: We solve the vacuum Einstein 
field equations in a 5-dimensional (5D) asymptotically 
flat spacetimc with an SO(3) symmetry. Since pertur- 
bations of 5D black strings violating this symmetry are 
stable and decay [l[ , we do not expect imposing this sym- 
metry qualitatively affects the results presented here. 

We use the generalized harmonic formulation of the 
field equations [31, and adopt a Cartesian coordinate 
system related to spherical polar coordinates via = 
{t,x,y,z,iv) = {t,rcos(j)sm9,r sin(j)sm9,r cosO, z). The 
black string horizon has topology 5*^ x R; {9, cj)) are co- 
ordinates on the 2-sphere, and z (iJu) is the coordinate in 
the string direction, which we make periodic with length 
L. We impose a Cartesian Harmonic gauge condition, 
i.e. VaV^x* = 0, as empirically this seems to result in 
more stable numerical evolution compared to spherical 
harmonic coordinates. The SO(3) symmetry is enforced 
using the variant of the "cartoon" method 120| described 
in [l9[ , were we only evolve ay = z = slice of the space- 
time. We further add constraint damping [2l| . which in- 



troduces two parameters k and p; we use {k,p ^ 1, —0.5), 
where a non-zero p is essential to damp an unstable zero- 
wavelength mode arising in the z direction. 

We discretize the equations using 4th order finite dif- 
ference approximations, and integrate in time using 4th 
order Runge-Kutta. To resolve the small length scales 
that develop during evolution we use Berger and Oliger 
adaptive mesh refinement. Truncation error estimates 
are used to dynamically generate the mesh hierarchy, and 
we use a spatial and temporal refinement ratio of 2. 

At the outer boundary we impose Dirichlet conditions, 
with the metric set to that of the initial data. These con- 
ditions are not strictly physically correct at finite radius, 
though the outer boundary is placed sufficiently far that 
it is causally disconnected from the horizon for the time 
of the simulation. We use black hole excision on the inner 
surface; namely, we find the apparent horizon (AH) us- 
ing a flow method, and dynamically adjust this boundary 




50 100 150 200 t/M 
FIG. 1: (Normalized) apparent horizon area vs. time. 



(the excision surface) to be some distance within the AH. 
Due to the causal nature of spacetimc inside the AH, no 
boundary conditions are placed on the excision surface. 

We adopt initial data describing a perturbed black 
string of mass per unit length M and length L = 20M w 
lALc {Lc is the critical length above which all pertur- 
bations are unstable). This data was used in Q and we 
refer the reader to that work for further details. 

We evaluate the following curvature scalars on the AH: 



K = /i?^^/12, 5 = 27 (UJ^I-^ - 1) + 1 



(1) 



where / = RatcdR''^"^, J = RabcdR'"^''^ Ref''^ and Rah 
is the areal radius of the AH at the corresponding point 
(though note that / and J are usually defined in terms 
of the Weyl tensor, though here this equals the Riemann 
tensor as we are in vacuum). K and S have been scaled 
to evaluate to {6, 1} for the hyperspherical black hole and 
black string respectively. 

Results: The results described here are from sim- 
ulations where the computational domain is (r, z) G 
([0, 320M] X [0,20M]). The coarsest grid covering the 
entire domain has a resolution of {Nr,Nz) = (1025,9) 
points. For convergence studies we ran simulations 
with 3 values of the maximum estimated truncation er- 
ror r: ["low", "medium", "high"] resolution have r = 
[to, to/&, To/ 64] respectively. This leads to an initial hier- 
archy where the horizon of the black string is covered by 
4, 5 and 6 additional refined levels for the low to high res- 
olutions, respectively. Each simulation was stopped when 
the estimated computational resources required for con- 
tinued evolution was prohibitively high (which naturally 
occurred later in physical time for the lower resolutions); 
by then the hierarchies were as deep as 17 levels. 

Fig. [T] shows the integrated AH area A within z S 
[0, L] versus time. At the end of the lowest resolution 
run the total area is A = (1.369 ± 0.005)^o[23|, where 
Aq is the initial area; interestingly, this almost reaches 
the value of 1.374Ao that an exact 5D black hole of the 
same total mass would have. Fig. [5] shows snapshots of 
embedding diagrams of the AH, and Fig. [3] shows the 
curvature invariants ([T]) evaluated on the AH at the last 
time step, both from the medium resolution run. 

The shape of the AH, and that the invariants are tend- 
ing to the limits associated with pure black strings or 
black holes at corresponding locations on the AH, sug- 
gests it is reasonable to describe the local geometry as 
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FIG. 2: Embedding diagram of the apparent horizon at sev- 
eral instances in the evolution of the perturbed black string, 
from the medium resolution run. R is areal radius, and the 
embedding coordinate Z is defined so that the proper length 
of the horizon in the space-time z direction (for a fixed 0, (^) 
is exactly equal to the Euclidean length of R(Z) in the above 
figure. For visual aid copies of the diagrams reflected about 
i? = have also been drawn in. The light (dark) lines denote 
the first (last) time from the time-segment depicted in the 
corresponding panel. The computational domain is periodic 
in z with period iz = 20Af ; at the initial (final) time of the 
simulation &Z = 20Af {&Z = 27.2Af). 
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FIG. 3: Curvature invariants evaluated on the apparent hori- 
zon at the last time of the simulation depicted in Fig. [2l The 
invariant K evaluates to 1 for an exact black string, and 6 for 
an exact spherical black hole; similarly for S ([1}. 



being similar to a sequence of black holes connected by 
black strings. This also strongly suggests that satellite 
formation will continue self-similarly, as each string seg- 
ment resembles a uniform black string that is sufficiently 
long to be unstable. Even if at some point in the cas- 
cade shorter segments were to form, this would not be a 
stable configuration as generically the satellites will have 
some non-zero z-velocity, causing adjacent satellites to 
merge and effectively lengthening the connecting string 
segments. With this interpretation, we summarize key 
features of the AH dynamics in Table HI 

We estimate when this self-similar cascade will end. 



Gen. 


U/M 


Rs,i/M 


Ls,i 1 Rs,i 


Us 


Rh.f/M 


1 


118.1 ± 0.5 


2.00 


10.0 


1 


4.09 ± 0.5% 


2 


203.1 ±0.5 


0.148 ±1% 


105 ± 1% 


1 


0.63 ± 2% 


3 


223 ± 2 


0.05 ± 20% 


« 10^ 


> 1 


0.1 - 0.2 


4 


« 227 


« 0.02 


^ 10^ 


> 1(?) 


? 



TABLE I: Properties of the evolving black string appar- 
ent horizon, interpreted as proceeding through several self- 
similar generations, where each local string segment tem- 
porarily reaches a near-steady state before the onset of the 
next GL instability, ti is the time when the instability has 
grown to where the nascent spherical region reaches an areal 
radius 1.5 times the surrounding string-segment radius Ra.i, 
which has an estimated proper length Ls,i (the critical L/R 
is ~ 7.2 Pi), ris is the number of satellites that form per 
segment, that each attain a radius Rhj measured at the end 
of the simulation. Errors, where appropriate, come from con- 
vergence tests. After the second generation the number and 
distribution of satellites that form depend sensitively on grid 
parameters, and perhaps the only "convergent" result we have 
then is at roughly t = 223 a third generation does develop. We 
surmise the reason for this is the long parent string segments 
could have multiple unstable modes with similar growth rates, 
and which is first excited is significantly affected by truncation 
error. We have only had the resources to run the lowest res- 
olution simulation for sufficiently long to see the onset of the 
4th generation, hence the lack of error estimates and presence 
of question marks in the corresponding row. 
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FIG. 4: Logarithm of the areal radius vs. logarithm of time 
for select points on the apparent horizon from the simulation 
depicted in Fig. (2] We have shifted the time axis assuming 
self-similar behavior; the putative naked singularity forms at 
asymptotic time t/M « 231. The coordinates at z = 15,5 
and 4.06 correspond to the maxima of the areal radii of the 
first and second generation satellites, and one of the third 
generation satellites at the time the simulation stopped. The 
value 2 = 6.5 is a representative slice in the middle of a piece of 
the horizon that remains string-like throughout the evolution. 



The time when the first satellite appears is controlled by 
the perturbation imparted by the initial data; here that 
is Tq/M « 118. Subsequent timescales should approxi- 
mately represent the generic development of the instabil- 
ity. The time for the first instability after that sourced 
by the initial data is Ti/M « 80. Beyond that, with 
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the caveats that we have a small number of points and 
poor control over errors at late times, each subsequent 
instability unfolds on a time-scale X k 1/A times that 
of the preceding one. This is to be expected if, as for 
the exact black string, the time scale is proportional to 
the string radius. The time <o of the end-state is then 
io « To + YZoTiX' = To+ ri/(l - X). For the 
data here, t^/M « 231; then the local string segments 
reach zero radius, and the curvature visible to exterior 
observers diverges. Fig. |4] shows a few points on the AH, 
scaled assuming this behavior. In the Rayleigh-Platcau 
analogue, the shrinking neck of a fluid stream has a self- 
similar scaling solution that satisfies r oc (Iq — t), or, 
dhir/d(— ln(to — t)) = —1, where r is the stream radius 
(see [i^l, and Q for extensions to higher dimensions); 
to within 10 — 20% this is the average slope we see (e.g. 
Fig. |4|) at string segments of the AH at late times. 
Conclusions: We have studied the dynamics of a per- 
turbed, unstable 5D black string. The horizon behaves 
similarly to the surface of a stream of low viscosity fluid 
subject to the Rayleigh-Plateau instability. Multiple gen- 
erations of spherical satellites, connected by ever thin- 
ner string segments, form. Curvature invariants on the 
horizon suggest this is a self-similar process, where at 
each stage the local string/spherical segments resemble 
the corresponding exact solutions. Furthermore, the time 
scale for the formation of the next generation is propor- 
tional to the local string radius, implying the cascade will 
terminate in finite asymptotic time. Since local curvature 
scalars grow with inverse powers of the string radius, this 
end-state will thus be a naked, curvature singularity. If 
quantum gravity resolves these singularities, a series of 
spherical black holes will emerge. However, small mo- 
mentum perturbations in the extra dimension would in- 
duce the merger of these black holes, thus for a compact 
extra dimension the end state of the GL instability will 
be a single black hole with spherical topology. 

The kind of singularity reached here via a self-similar 
process is akin to that formed in critical gravitational col- 
lapse f2^ : however, here no fine-tuning is required. Thus, 
5 (and presumably higher) dimensional Einstein gravity 
allows solutions that generically violate cosmic censor- 
ship. Angular momentum will likely not alter_this con- 
clusion, since as argued in [ij" 
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tion does not suppress the unstable modes, and moreover 
induces super-radiant and gyrating instabilities (26j . 
Acknowledgments: We thank V. Cardoso, M. Chop- 
tuik, R. Emparan, D. Garfinkle, K. Lake, S. Gubser, G. 
Horowitz, D. Marolf, R. Myers, W. Unruh and R. Wald 
for stimulating discussions. This work was supported by 
NSERC (LL), CIFAR (LL), the Alfred P. Sloan Foun- 
dation (FP), and NSF grant PHY-0745779 (FP). Simu- 



[1 
[2 

[3: 

[4: 

[s: 
[6: 

[7 

[9: 

[lo; 

[11 

[12 

[13 

[14 

[15 

[16 

[17 
[18 

[19 
[20 
[21 

[22' 
[23 

[24 
[25 

[26] 



R. Gregory and R. Laflamme, Phys. Rev. Lett. 70, 2837 
(1993). 

R. Emparan and H. S. Reall, Living Rev. Rel. 11, 6 
(2008); T. Earmark, V. Niarchos, and N. A. Obers, Class. 
Quant. Grav. 24, Rl (2007); B. Kol, Phys. Rept. 422, 
119 (2006). 

S. W. Hawking and G. F. R. Ellis, The Large Scale Struc- 
ture of Space- Time (Cambridge University Press, 1975). 
G. T. Horowitz and K. Maeda, Phys. Rev. Lett. 87, 
131301 (2001). 

T. Wiseman, Class. Quant. Grav. 20, 1137 (2003). 
S. S. Gubser, Class. Quant. Grav. 19, 4825 (2002). 

E. Sorkin, Phys. Rev. Lett. 93, 031601 (2004). 
M. W. Choptuik et al., Phys. Rev. D68, 044001 (2003). 
D. Garfinkle, L. Lehner, and F. Pretorius, Phys. Rev. 
D71, 064009 (2005). 

D. Marolf, Phys. Rev. D71, 127504 (2005). 
K. S. . Thorne, R. H. . Price and D. A. . Macdonald, 
"Black holes: The Membrane Paradigm," New Haven, 
USA: Yale Umv. Pr. (1986) 367p 

S. Bhattacharyya, V. E. Hubeny, S. Minwalla, and 
M. Rangamani, JHEP 02, 045 (2008). 
R. Emparan, T. Harmark, V. Niarchos, and N. A. Obers, 
JHEP 03, 063 (2010). 

V. Cardoso and O. J. C. Dias, Phys. Rev. Lett. 96, 
181601 (2006). 

V. Cardoso and L. Gualtieri, Class. Quant. Grav. 23, 
7151 (2006). 

J. Camps, R. Emparan and N. Haddad, JHEP 1005, 042 
(2010) 

J. Eggers, Rev. Mod. Phys. 69, 865 (1997). 
P. Kovtun, D. T. Son, and A. O. Starinets, Phys. Rev. 
Lett. 94, 111601 (2005). 

F. Pretorius, Class. Quant. Grav. 22, 425 (2005). 
M. Alcubierre et al.. Int. J. Mod. Phys. DIO, 273 (2001). 

C. Gundlach, J. M. Martin-Garcia, G. Calabrese, and 
I. Hinder, Class. Quant. Grav. 22, 3767 (2005). 
J. Eggers, Phys.Rev.Lett 71, 3458 (1993). 
U. Miyamoto, larXiv: 1007.4302 ^ [hep-th] : U. Miyamoto, to 
be published 

M. W. Choptuik, Phys. Rev. Lett. 70, 9 (1993). 
O. J. C. Dias, P. Figueras, R. Monteiro, H. S. Reall and 
J. E. Santos, JHEP 1005, 076 (2010); 

D. Marolf and B. C. Palmer, Phys. Rev. D 70, 084045 
(2004); V. Cardoso and J. P. S. Lemos, Phys. Lett. B621, 
219 (2005). V. Cardoso and S. Yoshida, JHEP 0507, 009 
(2005) 

[27] The error in the area was estimated from convergence at 
the latest time data was available from all simulations 



